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1 Introduction 

Modern semiconductor technology has in recent years made it possible to fabri- 
cate ultrasmall structures that confine electrons on scales comparable to their de 
Broglie wavelength. If the confinement is only in one spatial direction such systems 
are called quantum wells. In quantum wires the electrons can move freely in one 
dimension but are restricted in the other two. Structures that restrict the motion 
of the electrons in all directions are called quan tum dots. The number of electrons, 
N, in a quantum dot can range from zero to several thousand. The confinement 
length scales R±, i?2, can be different in the three spatial dimensions, but typ- 
ically i?3 <C R\ w i?2 ~ 100 nm. In models of such dots R3 is often taken to be 
strictly zero and the confinement in the other two dimensions is described by a 
potential V with V(x) — > 00 for |x| — > 00, x = (x l ,x 2 ) S R 2 . A parabolic poten- 
tial, V — ^uj\x\ 2 , is often used as a realistic and at the same time computationally 
convenient approximation. 

Quantum dots have potential applications in microelectronics and have been 
extensively studied both experimentally and theoretically. Apart from possible 
practical uses they are of great interest for basic quantum physics. Their parame- 
ters (strength and shape of the confining potential, magnetic field strength, number 
of electrons) can be varied in a controlled way and their properties can be studied 
by clever experimental techniques. This offers many possibilities to confront theo- 
retical predictions with experimental findings. There exist by now many excellent 
reviews on the physics of quantum dots, e.g. [jjj-j^]. In the present contribution 
the focus will be on some theoretical aspects that are only partly covered by these 
reviews, in particular on rigorous limit theorems ]To| which apply to dots in 
high magnetic fields and/or with high electron density. 

A quantum dot with N electrons is usually modeled by a Hamiltonian of the 
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following form, acting on the Hilbert space 

N 

H N = f\L 2 (R 2 ;C 2 ) (1) 
1 

appropriate for two dimensional Fermions of spin 1/2: 

JV 

H N = J2H { / > + Yl W( Xi - Xj ), (2) 

j=l l<i<j<N 

where x t <E R 2 , i = 1, . . . , N and Hi = l<8>---<8).ffi®---®l (Hi in the j-th 
place) with the one-body hamiltonian 

h 2 / e \ 2 / Tie \ 

ffi = 2^: ( iV ~ ^ AW ) + v & + 5 * ( 2^ J - (3) 

Here = j(—Bx 2 1 Bx 1 ) is the vector potential of a homogeneous magnetic 

field of strength B in the x 3 -direction, V is the confining potential, assumed to 
be continuous with V(x) — > oo for \x\ — > oo, and S3 is the spin operator in in - 
direction. The parameters m* and are respectively the effective mass and the 
effective g-factor of the electrons, while m e and e are the bare values of the electron 
mass and electric charge, and K and c have their usual meanings. The constant 



\2m c J \m* 2 

has been introduced in (3) for convenience: Subtraction of CB has the effect that 
the spectrum of the kinetic energy operator (including spin) iJ}" 11 = Hi — V starts 
at zero for all B, even if m* ^ m c and <?» 7^ 2. 

The interaction potential W represents the Coulomb repulsion between the 
electrons, modified by the properties of the surrounding medium. Usually it is 
simply taken to be 

W(xi - xj) = el\xi - Xj\~ x (5) 

where e* = e/^/e with e the dielectric constant, but some regularization of the 
bare Coulomb potential, e.g., 

W{ Xl - Xj ) = el [(\ Xi - Xj \ 2 + dl)- 1 ' 2 - (\xi - x,\ 2 + S 2 _ )- 1 / 2 ] (6) 

with5_ > 6+ > @, or even a potential that depends not just on the differences 
Xi — Xj, may sometimes fit the effective interaction better. For the proof of some 
of the theorems below the important property of (||) is that W is repulsive, of 
positive type and tends to zero at infinity; these features are shared by (||). 

Writing the Hamiltonian in the above form is, of course, an approximation, 
because the effect of the medium on the electrons is only taken into account 



2 



through the modification of the parameters from their bare values. The size of 
a quantum dot (~ 100 nm) is however usually much larger than the lattice con- 
stant of the medium where it resides (< 0.5 nm), so this approximation is usually 
a good one. 

Quantum dots, especially such with few electrons, are sometimes referred to 
as artificial atoms with V playing the role of the attractive nuclear potential in real 
atoms. The analogy is not perfect, however, because V is regular around the origin 
in contrast to the potential from an atomic nucleus, and also because the electron 
interaction is the three dimensional Coulomb potential (|5|) (or modified Coulomb 
potential @), while the motion is (essentially) restricted to two dimensions. But 
in many respects quantum dots can indeed be regarded as artificial atoms, with 
an important additional aspect: The effective parameters are to a certain extent 
tunable and may differ appreciably form their counterparts in real atoms. 

In a quantum dot the natural atomic unit of length is a* = eh 2 /(m*e 2 ). 
Compared with the usual Bohr radius, ao = h 2 / (m c e 2 ) = 0.53 x 10 _1 nm, the 
length a* is typically large, e.g., a* « 185 ao ~ 10 nm in GaAs. The natural 
energy unit is E* = e 2 /a* — e^m^/Ti 2 , and in GaAs E* ~ 12 meV, which should 
be compared with E — e 2 /a = e A m lc /fi 2 — 27.2 eV, i.e., J?,»4x \Q~ i E . 

The natural unit, £?*, for magnetic field strength is the field at which the 
magnetic length Ib = he/(B 1 ^ 2 c) equals a*, or equivalently, at which the cyclotron 
energy TieB^/m^c equals E*. Hence = (a Q /a*) 2 B , where B — e 3 m 2 c/h 3 = 
2.35 x 10 5 T is the value corresponding to free electrons. If ao/a* is small, B* can 
be much smaller than Bo. Thus B* » 7 T in GaAs. This opens the very interesting 
possibility to study in the laboratory magnetic effects, whose analog for real atoms 
require field strengths prevailing only on neutron stars. 

On the experimental side the main techniques for studying quantum dots 
are charge transport and capacitance spectroscopy ||, |[3|) and optical far 



infrared spectroscopy (p4|, 15 ). The former is in particular suited for measuring 
the N dependence of ground state energies, but excited states can be investigated 
as well. The applicability of optical spectroscopy is to a certain extent limited 
by Kohn's theorem, to be discussed below, but refined techniques allow also to 
infer many properties by this method. Altogether it is fair to say that the en- 
ergy spectrum of quantum dots and even some aspects of the corresponding wave 
functions are experimentally accessible with considerable precision. See also [[l6| 
for recently discovered effects in charge transport spectroscopy that wait for an 
adequate explanation. 

On the theoretical side the spectral properties of the Hamiltonian (^) have 
been studied by a variety of methods, which can be roughly divided into the 
following categories: 

• Exact analytic solutions 

• Rigorous limit theorems 

• Numerical diagonalizations 
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• Hartree and Hartree-Fock approximations 

• Variational calculations 

• Density functional methods 

In the condensed matter literature the last four categories are by far the most 
prominent, but the present survey is mostly concerned with the first two, which lie 
within the realm of mathematical physics. It is impossible here to do any justice to 
the extensive physics literature on the theory of quantum dots, but an annotated 
bibliography of some representative references will be given in the last section. 

2 Exact solutions 

2.1 The Fock-Darwin spectrum 

From now on units arc chosen such that Ti = e* = to* = = 1. The Hamiltonian 
(||) can then be written 

Hi = \ (iV - A{x)f + V(x) + -yS 3 B - |(1 - h\)B (7) 

with 7 = <7„<TO*/(2TOe). For a confining potential of the harmonic oscillator form 

V(x) = ± W |x| 2 , (8) 

the eigenvalues and eigenfunctions of 

H? b = ±{iV-A(x)) 2 + ^ 2 \x\ 2 (9) 

were determined by Fock already in 1928 [[l7|, and also by Darwin in 1930 Ell . 
The Fock-Darwin spectrum of (0) consists of the eigenvalues 

e£? = (2fc + |l| + l)fi-§IB (10) 

with 

f7=(IS 2 +c 2 ) 1/2 , (11) 
k = 0, 1, 2 . . ., / = 0, ±1, ±2, . . .. The corresponding eigenfunctions are 

^(i) = (const.) exp^)?- 1 ' 1 exp(-f7r 2 /2)L^ {Qr 2 ) (12) 

where x = (r cos (p,r sin<p) and is an associated Laguerre polynomial. These 
functions are also eigenfunctions of the angular momentum L3 in the a; 3 -direction 
with eigenvalue I. Eigenfunctions with the same value of 

n=k + \{\l\-l) (13) 
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are grouped together in a Fock-Darwin level (FDL). In the limit lu/B — > the 
eigenvalues in a FDL coalesce and a FDL becomes identical to a Landau level 
(LL) with the eigenvalues 

£ £ = (n+i)5 (14) 

and eigenfunctions 

ipl tl (x) = (const.) exp(iZ^)r |(| exp(-Br 2 )Lj £ ' 1 (2Br 2 ) (15) 

The degeneracy of a LL per unit area is Bj (2n). 

If the interaction W between the electrons is ignored the FD spectrum to- 
gether with the Pauli principle, taking spin into account, completely solves the 
eigenvalue problem for Hjy in the case of the quadratic potential (^) . This approx- 
imation even fits some experimental data quite well fl , [0 . When the interaction 
W is taken into account this picture has to be modified, of course. For a quadratic 
confining potential, however, the FD spectrum continues to apply to the motion 
of the center of mass, independently of the interaction. This simple, but impor- 
tant fact jni, |Q goes under the heading Kohn's theorem. The proof is essentially 
contained in the identity 

N I N \ 2 

N J2 x '= 5>i +E^-^-) 2 ' ( 16 ) 

3=1 \3=1 / i<J 

for it implies that Hm can be written 

H N =H™+H s v' m + H^ (17) 

where Hff operates only on the relative coordinates Xi — Xj , while all dependence 
on the center of mass coordinate X = (x\ + • • • + Xn)/N is contained in 

H™ = i (iV x - NA{X)f + f W 2 |X| 2 . (18) 

The spectrum of 7?^ M is exactly the same as that of H ° rb , independent of TV and 
14^. In the dipole approximation the radiation field couples only to X and in this 
approximation optical FIR spectroscopy thus probes only the FD spectrum. 

There is a further instance where the FD eigenfunctions play a role even 
for W ^ 0. Let IIo D denote the projector on the subspace oi Hn generated by 
eigenfunctions with FD index n = and complete polarization, i.e., spin magnetic 
moment in the direction of the field. Consider the operator 

H' N = n^H N U F B . (19) 

For large B the spectral properties of this operator can be expected to approximate 
those of Hm- If W — the ground state of this operator has the orbital wave 
function 

lT ,mdd / FD A / FD A A / FD / or ,\ 
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This eigenfunction is called the maximum density droplet [El] because the electrons 
are as "tightly packed" as possible around the origin. It is an amusing observation 
that \I/ mdd remains an exact eigenfunction of H' N also for W 7^ 0. The proof is very 
simple: H' N and the angular momentum operator L3 on 7ijv commute and have 
discrete spectrum. Hence for every eigenvalue of L3 there is a corresponding eigen- 
function that is simultaneously an eigenfunction of H' N . But the lowest eigenvalue, 
N(N— 1)/2, of L3 in the subspace of Hat generated by the FD functions with n = 
is nondegenerate, and the corresponding eigenfunction v|/ mdd must hence also be 
an eigenfunction of H' N . There is numerical ^2), [§3|) and even experimental 

p3 evidence that \I/ mdd is a ground state of H' N for some range of values of B. 

2.2 Analytic solutions for N > 1 

For N = 2 the Hamiltonian for the relative coordinate x = x± — X2 is 

#2 Cl = (iV - \A{x)f + \u 2 \x\ 2 + W(aO (21) 

In the case of a pure Coulomb interaction, W(x) = l/\x\, explicit formulas for 
eigenfunctions and eigenvalues of (|T]) have been found by Taut ■ His approach 
is based on an ansatz for the wave functions of the form 

ip(x) = exp(im^) exp(-p 2 / 2 )/° |m| - p (p) ( 22 ) 

where P is a polynomial in p = (Q/2) 1 / 2 r, with £1 = (jB 2 + a; 2 ) 1 ^ 2 , <p is the 
angular variable and m € Z. Solutions of the form ( p2| ) with P a polynomial do not 
exist for arbitrary values of B, but Taut's method produces at least eigenfunctions 
and eigenvalues for a countable infinity of values of which accumulate at 0. In 
fact, an ansatz like ( p2| ) in the eigenvalue equation i?| e V = Eip with P a power 
series, 

P{p) = Y / ^p\ (23) 

leads to 

a u = F{\m\,u,E')a (24) 
with a certain recursively computable function F and where E' is related to E by 

£ = ifiS' - \mB. (25) 

The condition that a v vanishes for all v larger than some positive integer n is 
equivalent to the two equations 

F(\m\,n,E',n) = (26) 

and 

E' = 2(\m\+n). (27) 
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For given n and m this gives one or more acceptable values for Q and corresponding 
energy values 

E = \{n + \m\)n-\mB. (28) 

The solutions found in this way are not necessarily ground states of (^l|) but the 
position of E in the spectrum can be inferred from the number of nodes of the 
corresponding wave function. 

The solutions of Taut seem so far to be the only known exact solutions for 
N = 2 and the Coulomb interaction (||). They are limited to the special values of 
VL defined by (^6|) and ( p7| ) . For W of the inverse square form 

W{x)=a\x\- 2 (29) 



on the other hand, the Hamiltonian ( |2l| ) can be exactly diagonalized fo all B p6| . 
In fact, addition of ( p9|) merely modifies the centrifugal term in the radial part of 
the FD Hamiltonian (g) and we obtain as eigenvalues of ( ^lf ) 

E= [2n + /x+l]Q- \mB (30) 

with \x = (m 2 + a) 1 / 2 not necessarily an integer, and the eigenfunctions 

r/,{x) = exp(im^) exp( V /2)p"L£(p 2 ). (31) 



It should be noted that the inverse square form (M) for the effective interac- 
tion between the electrons is not necessarily less realistic than the pure Coulomb 
repulsion (|J). In fact, the form (|^) of the interaction, that is motivated by the situ- 



ation in real dots 11 , has an inverse square decrease for large separation. At small 
separation, on the other hand, the effective interaction may be less singular than 
(p|) . It is therefore not entirely academic to consider also a harmonic interaction 
p7fl of the form 

W(x i -x j ) = 2W o -±0\x i -x j \ 2 (32) 

with positive parameters Wo and f3. For this case one can even solve the problem 
for all N exactly, using ([If]) to decouple the oscillators. The result for the ground 
state energy E Q (N, B) of H N with N > 2 and 7 = (for simplicity) is |27j 

E Q (N, B) = fi + \(N - 1)(N - 2)Q (N) - \N(N + 1)5 + N(N - 1)W (33) 

with Q (N) = (n 2 - NP 2 ) 1 ' 2 . It is assumed that uj > iV 1 / 2 /?, so fl (N) € R + for 
all B. The corresponding wave function for the relative motion is 

^j(x) = Y[[ Zij exp(-r!o|%-| 2 /2iV)] (34) 

i<j 

where — {x\ — xj) — i(xf — x 2 ) are the relative coordinates regarded as points in 
C. The form ( |3^ ) of the repulsion, of course, quite wrong for large separation, but 
this error is to some extent counterbalanced by the confining potential. For a com- 
parison of the solutions for different 1-F's and their confrontation with numerical 
calculations and experiments we refer to B. 
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3 Rigorous limit theorems 



While exact solutions of the eigenvalue problem for Hn with W of the form or 
are not available for N > 2, it is possible to analyze at least the ground state 
properties for large large B or large N by minimizing simple functionals of the 
electron density. In this analysis, which implies a drastic reduction of the quantum 
mechanical A^-body problem, it is not necessary that the confining potential V 
has the quadratic form (^) . From now on V will stand for an arbitrary continuous 
function of x G R 2 tending to oo for |x| — > oo; when additional properties are 
required these will be explicitly stated. It is no restriction to assume that V > 0. 
To be able to consider variations of the strength of the potential at fixed shape, 
we write 

V{x) = Kv(x) (35) 

with a coupling constant K. The interaction W will for definiteness be assumed 
to be pure Coulomb with e* = 1, i.e., 

W(x-y) = \x-y\- 1 , (36) 

but other repulsive potentials of positive type could be treated similarly. The 
quantum mechanical ground state energy is 

E®(N,B,K) = {$ ,H N y), (37) 

where is a normalized ground state of Hjy. The corresponding ground state 
electron density is 

Pn,b,k( x ) = J a ii x 2,o-2] ■ ■ ■ ,x N ,a N )\ 2 dx 2 ■ ■ -dx N . (38) 

spins (7k— ±1/2 

We are concerned with the asymptotics of these quantities when one or more of 
the parameters N, K and B tends to oo with v fixed. 

The large B limit at fixed N and K is easiest and will be considered first. 

3.1 High field limit 

In the lowest Landau level (i.e., for n = 0) the wave functions ( fllf ) are essentially 
localized on scale B^ 1 / 2 , and the quantum mechanical kinetic energy vanishes 
after the spin contribution and the subtraction of ^(1 — I7I) in (]?]) have been taken 
into account. In the limit B — > 00 it can therefore be expected that a classical 
model of N point particles with the energy function 

N 

£ F [xi,...,x N ] =^2v(x i ) + ^2\x i -x j \- 1 (39) 

i—1 i<j 
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describes the ground state energy correctly. This is indeed the case: Defining 

E P (N,K) =M£ p [ Xl ,...,x N ] (40) 



we have 

Theorem 3.1 (High field limit.) For N and K fixed, 

lim E Q (N, B, K) = E P (N, K). (41) 

B^oo 

Proof: The lower bound for in the proof of this theorem is trivial: Since 
the kinetic energy part of (@) is > 0, E p < E®. The upper bound is a simple 
variational calculation: Since V(x) tends to oo for \x\ — > oo and \xi — Xj\ — > oo 
for \xi — Xj \ — > 0, the infimum of £ p is obtained at some point (x\, . . . , xn) € H 2N 
with Xi Xj for i ^ j. We test the Hamiltonian with the wave function 

^(xi . . . , xn) = exp(ixi x xi/ B)tpQ (xi— xx) A. . .Aexp(ixN x xn/B)iPq (xn — x^) 

(42) 

A indicates antisymmetrization in the coordinates (xi, . . . , Xjsr), and the vector 
product x of two points in R 2 is regarded as R valued. The spins are totally 
aligned and are therefore omitted in the notation. The kinetic energy of W is 
zero, because the wave functions exp(kc x Xj/B)^ (x — Xj) belong to the lowest 
Landau level. The potential energy part approximates £ p [xi, . . . , xn} arbitrarily 
well as B — > oo, because \ipQ Q (x — Xj)\ 2 is essentially localized within a radius 
~ B^ 1 / 2 around Xj and tends to a delta function as B — > oo. 

3.2 Large N limits 

Before stating the limit theorems formally let us briefly discuss their heuristic 
basis. 

For B — each of the N electrons in a dot of radius R occupies a "private 
room" of spatial extension R/N 1 / 2 , because of the Pauli principle. Hence the 
kinetic energy of an electron is £kin ~ N/R 2 . The potential energy due to the 
confining potential is e CO nf ~ Kv{R) and the repulsive energy due to the other 
electrons £ rcp ~ N/R. The radius R can be estimated by minimizing the sum £kin+ 
£conf + £rep- If K/N is kept fixed while N — > oo, all these terms are proportional 
to N. The radius is therefore independent of N an the mean distance between the 
electrons is ~ RN -1 / 2 . In this limit the density tends to infinity. The ground state 
energy of the dot, i.e., the total energy E® of the N electrons, is ~ TV 2 , and this 
is large compared to the exchange/correlation energy ~ N ■ (TV -1 / 2 ) -1 = iV 3 / 2 . 
One may therefore expect that a Thomas-Fermi theory pj| captures the leading 
asymptotics for N, K — > oo with K/N fixed. 

A magnetic field B will not influence the asymptotics as long as B <C N, 
because the energy differences between Landau levels are much smaller than the 
other energy contributions. If B ~ N, on the other hand, all Landau levels have 
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to be taken into account, and the Thomas- Fermi theory has to be modified (fLl[, 
p9[ ) . For B ^> N the electrons will essentially all sit in the lowest Landau level and 
the kinetic energy contribution vanishes. The radius is essentially independent of 
B and N (as long as B ^> N) and slightly smaller than in the B = case, because 
the positive kinetic energy term £ki n ~ N/R 2 is now missing. The ground state 
energy is ~ iV 2 as before. Theorem 3.1 applies, but the large N limit leads to 
an additional simplification: The point charges can be replaced by a continuous 
distribution because the interelectronic distance ~ RN -1 / 2 tends to zero. 
The asymptotics is described by the energy functional of a charged fluid in the 
confining potential V. 

For a homogeneous potential V of degree s > 1, i.e., V(Xx) — X s V(x), this 
last mentioned "classical" energy functional describes also the large N limit for 
arbitrary B if K/N — > 0, in particular if K is fixed while N — > oo. The radius 
R ~ (N/K) l ^ s+1 ^ tends in this case to oo. The density may go to zero or to infinity. 
The kinetic energy per particle ~ N/R 2 ~ N(K/N) 2 ^ S+1 " > is in both cases small 
compared to the other contributions ~ KR S ~ K(N/K) S ' Smearing out the 

point charges brings in an error ~ N 1 / 2 /R per particle, but it is small compared 
to the total Coulomb repulsion ~ N/R. 

These heuristic considerations will now be turned into precise statements. We 
define three functionals of the electron density p £ Li(R 2 , dx), p > as follows. 

The 2D Thomas-Fermi functional is 

£ TF [p;if] = 7T J p(x) 2 dx + J V(x)p(x)dx + D(p,p) (43) 

with D(p, p) = | J J p{x)\x — y\~ 1 p(y)dxdy. All integrals are over R 2 . The np 2 in 
the first term is just the kinetic energy density of a two dimensional noninteracting 
electron gas of density p at B = 0. 

The 2D magnetic Thomas-Fermi functional is defined as 
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[p; B,K} = J ]B (p(x))dx + J V(x)p(x)dx + D(p, p) (44) 



where js is a piece wise linear function representing the kinetic energy density of 
the electrons in a magnetic field B, taking all Landau levels (including spin) into 
account. If 7 in ((7J) were 0, the derivative j' B = djs/dp would just be the step 
function 

j' B (p)=B[2irp/B] (45) 

where [t] denotes the integer part of a real number t. This, together with js(0) = 
fixes js(p). As B — > 0, - * np 2 . If 7 7^ 0, then j B is a slightly more compli- 
cated step function, cf. Eq. (2.3) in |J, but since the explicit form is not significant 
we refrain form stating it here. When a concrete j'b is needed for discussion pur- 
poses we shall stick to the simplest case (^5|). 
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Finally, the classical energy functional is defined as 

£ c [p;K} = f V(x)p(x)dx + D(p,p). (46) 



All three functionals are convex, and for each there is a unique nonnegative 
density that minimizes the functional under the constraint J p(x)dx — N; this is 
discussed in ||. 

We denotes the minimizing densities and the corresponding energies by pJf K 
and E TF (N,K) for the TF functional @, by Pn T b f k and E MTF (N,B,K) for 
the magnetic TF functional Q) , and by p^ K and E c (N, K) for the classical 
functional ( (f6| ) . In the last case some regularity of V is needed in order to ensure 
that Px K is a function and not just a positive measure; a sufficient conditions is 
that V satisfies locally an estimate of the form 

| W(x) - VV(y)\ < (const.) \x - y\ a (47) 

with some a > 0. 

The TF functional J43| ) and the classical functional (^6|) are both limiting 
cases of the MTF functional (]44|) , for B — > and — > oo respectively. More precisely, 

\im o E MTF (N,B,K) = E TF (N,K) (48) 

limpTJU = ^ (49) 

and 

lim E MTF {N,B,K) = E C (N,K) (50) 

hm p™ TF K = (51) 

B—>oo 

The limit for the densities should be understood in the weak L\ sense, but for 
special V much stronger convergence may hold. For instance, if V is monotonically 
increasing with |x|, p% K is a bounded function, and p^ K — P^ T b k ^ ot sufficiently 
large B, because Jb{Pn k) — ior B > 27r|| j o c || 00 . 

The MTF theory has two nontrivial parameters because of the scaling re- 
lations 

E MTF (N,B,K) = N 2 E MTF (1,B/N,K/N) (52) 
= Npf™N,K /N {x)- (53) 

Corresponding relations (without B) hold for the TF theory and the classical 
theory, and also for E p . 

A further important property of the densities is their compact support: For 
fixed K/N the minimizers of £ TF , £ MTF , £ c and also of £ p have support in a disc 
whose radius is uniformly bounded in TV and B (Lemma A.l in [§)). 
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Each minimizcr satisfies a variational equation, which in the case of the MTF 
theory has an unusual form, since it consists really of inequalities. To state it 
compactly it is convenient to modify the definition ( f45| ) slightly and regard j' B as 
an interval valued function if 2np/ B is an integer, namely, if 2np/ B = n, then 
j' B is the closed interval [(n — l)B,nB}. The MTF equation that is satisfied by 
pMTF can t nen kg wr itten 

,-vw-,.M- « (54) 

with a unique p — p(N, B, K ). Such generalized variational equations have been 
studied by Lieb and Loss [j30| . 

If the potential is quadratic, V(x) = K\x\ 2 , there is an explicit formula (|p8|, 
for the minimizer for £ c , which is equal to p MTF for B sufficiently large: 

c (x) _ f &N\y/l-\\x\* if \x\ < A- 1 , . 

Pn,k( x )-\ o if | x | > \-i i 55 ) 

with A = (8K/3ttN) 2 / 3 . The density profile has the shape of a half ellipsoid with 
a maximum at x — 0. Note the difference between the two dimensional case con- 
sidered here, and three dimensional electrostatics: In three dimensions the density 
would be homogeneously distributed in a ball. 

The criterion for p MTF = p G is that js(/O C (0)) = 0, which holds if 

B > (6/3 2 /V/ 3 )X 2 /3at1/3. (56) 

Numerically computed profiles of the minimizers p MTF and the corresponding 
effective potentials 

V eS (x)=V(x)+p MTF *\x\- 1 (57) 

with V(x) = K\x\ 2 are shown in Fig. 1. The computations were carried out by 
Kristinn Johnsen. 

At the highest value of the field (Fig. 1(a)) condition ( |56| ) is fulfilled and p MTF 
has the form (|55|). On the support of p MTF — p c we have V e g(x) =constant= p. 

When the field is gradually turned down the maximal density B/(2n) of 
electrons that can be accommodated in the lowest Landau level goes down also. 
Condition (|5^) no longer holds, i.e., the density p c near the center is higher than 
B/(2tt) and charges have to be moved into other states in phase space. If B is only 
slightly smaller than the value given by (|5(j) (Fig. 1(b)) it would cost too much 
energy to bring the electrons near the origin into the next Landau level and it pays 
to move them spatially away from the center, because the potential energy increase 
is less than B. Hence in a certain range of B values, the density near the center is 
locked at the value B/(2n) ( "incompressible" domain). The effective potential is 
no longer constant in this domain. In the complementary "compressible" domain, 
on the other hand, the density is below the critical value Bj (2tv), and tends to zero 
in such a way that the effective potential stays constant. Reducing the field strength 
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Figure 1: Density profiles and effective potentials for the MTF theory at different 
magnetic field strengths, calculated for N — 50 and V{x) = K\x\ 2 with K = 1,7 
meV and the material parameters of GaAs. (a) B = 8 T, (b) B = 7 T, (c) B = 2 
T, (d) B = T. 
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further brings more Landau levels into play (Fig. 1(c)). Incompressible domains, 
where the density is an integer multiple of B/(2w), alternate with compressible 
domains, where the effective potential has a constant value. When B — > the 
profile becomes indistinguishable from the smooth profile of p TF (Fig. 1(d)). It is 
interesting to note that the alternation of compressible and incompressible domains 
in moderate magnetic fields may account for some fine structure in the charge 
transport spectroscopy of quantum dots with a large number of electrons pi|. 

The basic limit theorem || that relates the energy functionals ([i"3"|)-(|46D to 
the quantum mechanical ground state of Hn is as follows: 

Theorem 3.2 (High density limit.) Let N — > oo with K/N fixed. Then, uni- 
formly in B/N, 

E Q (N,B 1 K)/E MTF (N,B,K)^l (58) 

and 

N'^bM*) "> /%R/n,k/nW (59) 
in weak L 1 sense. Moreover, if B/N — > 0, then 

E Q (N,B,K)/E TF (N,B,K) -> 1 (60) 

N^Pn^.k^) -> Pi F b/n,k/n( x ), ( 61 ) 

and if B/N — > oo, then 

E Q (N,B,K)/E C (N,B,K) -> 1 (62) 

n ^P%,b,k( x ) -> Px,b/n,k/n( x )- ( 63 ) 

According to this theorem there are thus three asymptotic regimes for quan- 
tum dots as N and K tend to oo with K/N fixed: B < N, B ~ N and N < B. 
This should be compared with the more complex situation for three dimensional 



natural atoms in strong magnetic field, where there are five regimes 32 , [B3| 
for N -> oo with Z/N fixed (Z = nuclear charge): B < TV 4 / 3 , B ~ iV 4 / 3 , 
iV 4/3 «B« iV 3 , B ~ TV 3 , TV 3 < B. 

For homogeneous potentials a stronger asymptotic result holds, for K/N may 
tend to zero as N — > oo. 

Theorem 3.3 (Homogeneous potentials.) Assume that V is homogeneous of 
degree s > 1, i.e., 

V(Xx) = X s V{x). (64) 

Then 

lim E Q {N,B,K)/E MTF (N,B,K) = 1 (65) 

uniformly in B and in K as long as K/N is bounded above. Moreover, if K/N — > 
as N — > oo, t/ien 

lim E Ql {N,B,K)lE c {N,K) = 1 (66) 

iV— >oc 

uniformly in B. 
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We shall now discuss briefly the main techniques used for the proof of these 
theorems. As usual it is sufficient to prove the limit theorems for the energy, 
because the corresponding results for the density can be obtained by variation 



with respect to the potential V. The basic result is thus Eq. (58); the other limit 
theorems follow by ([l8])-([)l]). One has to prove upper and lower bounds for the 
quantum mechanical energy E® in terms of the energy £ MTF , with controllable 
errors. 

The upper bound is obtained, using the variational principle of H] , by testing 
-Hjv with a suitable one particle density operator. Its kernel in the space and spin 
variables x, a has the form 

fC(x, a; x', a 1 ) = f u (u)Yl vu {x, a; x', a')d 2 u (67) 

where the sum is over all Landau levels and f v {u) is the filling factor of the z;-th 
Landau level at point u when the density is p MTF (u). The kernel H„ u (x, a; x', a') 
is obtained from the kernel Ii v (x, cr; x' , a 1 ) of the projector on the v-th. Landau 
level by localizing around u with a smooth function g of compact support, i.e., 

ll vu (x, cr; x', cr') = g(x - u)U v (x, cr; x', a')g(x' - u). (68) 

This operator is positive and approximately a projector, localizing simultaneously 
in space, i.e., around u, and in the Landau level index v. By letting the support of 
g shrink with N more slowly that the average electron spacing TV -1 / 2 , the error 
terms in the estimate above for — E MTF are of lower order than TV 2 , which is 
the order of E MTF . 

The lower bound for is proved separately for large B and for small B. For 
large B, i.e., B ^> N, one starts with the obvious estimate E® > E p . One then has 
to compare E p with E c , i.e., the energy of point charges with those of smeared 
charges. Since the electron distance is ~ N~ x / 2 the self energy of a smeared unit 
charge is ~ TV 1 / 2 . Hence an estimate 

E p (N, K) >E C (N,K)- bN 3/2 (69) 

with b depending only on K/N is to be expected, and this can indeed be proved, 
using an electrostatic lemma of Lieb and Yau |3q| . 

The lower bound for small B, i.e., B <C N or B ~ N, requires an estimate 
on the indirect Coulomb energy, that is derived in essentially the same way as a 
corresponding inequality in p8[ , using the positive definiteness of the Coulomb 
interaction §, cf. also |§. 

Lemma 3.1 (Exchange inequality in 2 dimensions.) 

E / i*i 2 E^-^i _1 ^^*"°*)- 192 ( 2 ^ 1/2 /^ /2 - ( 7 °) 

spins R2N i<j R2 
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In order to control negative term ~ J p 3 / 2 on the right side of ( f70| ) a lower bound 
on the kinetic energy is needed. This in turn is derived from a two dimensional 
magnetic Licb-Thirring inequality, which has to be proved in a slightly different 
way from the corresponding inequality in J33|, because there is no kinetic energy 
associated with a motion in the x 3 direction. The following inequality is adequate 
for the present purpose, but sharper Lieb-Thirring type inequalities that hold even 
for inhomogeneous fields have been derived by Erdos and Solovej |36]|, ]37j |. 

Theorem 3.4 (Lieb-Thirring inequality in 2 dimensions.) Let U be locally 
integrable and let ei(U), e2(U), . . . denote the negative eigenvalues (if any) of the 
Hamiltonian ±(iV - A) 2 + S 3 B - U. Define \U\ + {x) = %[\U(x)\ + U(x)]. For all 
< A < 1 we have the estimate 

EMEOI^A- 1 ^ / \U\ + (x)dx + l(l-X)- 2 [ \U\\{x)dx. (71) 

By a Legendre transformation with respect to U it follows from ( fn] ) that for 
all < A < 1 the kinetic energy of a state *f> is bounded below by |(1 — 
A) 2 — A _1 B/27r] 2 h . It is then possible, for B/N smaller than a certain critical 
value depending on v, to chose an N dependent e > in such a way that e — ► 
as N — > oo, but eT^p — (const) J p 3 / 2 > for all iV-particle states 



4 Other approaches 

The rigorous results presented above concern mainly (but not exclusively) the 
extreme cases of very few (N = 1 or N = 2) or very many (N — > oo) electrons. 
These cases play a similar role as the hydrogen atom and the Thomas-Fermi atom 
do in ordinary atomic physics, i.e., they set a standard that can be used as a 
starting point of various approximation schemes, or as a test for such schemes 
whose connection with the original Hamiltonian (||) may not be entirely clear. 
The physics literature on quantum dots is by now quite extensive, and many 
approaches have been used for gaining insight where rigorous results are not yet 
available. Here it is only possible to mention the main methods and give a list of 
some references that are representative for the approaches of condensed matter 
physicists to these problems and from it further sources can be traced. See also [|| 
for a more comprehensive list. 

For small N a direct numerical diagonalization of -ffy is possible and has 
been carried out, e.g., in @, @, j4§, @, @, for various values of N < 8. It 
is, of course, necessary to restrict to a finite dimensional subspace of the full 
Hilbert space, and the error made in this step is seldom estimated rigorously. By 
the mini-max principle, however, the computed values give at least upper bounds 
to the true eigenvalues. One of the features studied by this method is the orbital 
angular momentum and spin of the ground state, and analogs of Hund's rules from 
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atomic physics, as well as oscillations between triplet and singlet states for N = 2 
as B is varied have been seen in the calculations pd| . 

For N > 10 numerical diagonalization of the Hamiltonian is at present hardly 
feasible and resource is taken to other techniques of many body theory like Hartree 
and Hartree- Fock approximations (e.g., p2| , p3[ , [ p2[ , Esl), perturbation theory 
Q, variational methods (e.g. Q), quantum Monte Carlo methods (Q, Q), 
and (current) density functional theory (Q, p3[ , |5l| , Q). In quantum dots 
with a moderate electron number correlations play a much greater role than in 
natural atoms because the confining potential is usually quite shallow around the 
origin and the density may be low. (By contrast, one of the main steps in the 
proof of Theorem 3.2 is to show that exchange and correlation effects vanish in 
the high density limit for arbitrary magnetic fields.) Strong corelations together 
with strong dependence on magnetic fields make up much of the special flavor 
of quantum dot physics and it remains a challenge for mathematical physics to 
capture these effects in a rigorous way. 
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